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Abstract 

The  current  distribution  over  the  plate  surface  in  lead-acid  cells  in  the  course  of  discharge  was  determined  mathematically  by 
using  the  equivalent  circuit  method.  The  dependence  of  the  internal  cell  resistance  on  the  current  and  charge  passed  was  determined 
by  measurements  on  a  laboratory  cell.  Six  cell  variants  were  considered  differing  by  the  location  of  tabs  serving  as  current 
terminals.  The  results  are  presented  in  the  form  of  3D  diagrams  at  various  states  of  discharge.  To  make  the  current  distribution  nearly 
uniform,  extended  current  tabs  located  at  opposite  ends  of  the  plate  electrodes  were  proposed.  ©  2002  Elsevier  Science  B.V.  All  rights 
reserved. 
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1.  Introduction 

At  higher  discharge  rates,  the  performance  of  lead-acid 
batteries  may  be  adversely  affected  by  non-uniform  current 
distribution  over  the  plate  surface  owing  to  Ohmic  losses  in 
the  solid  phase.  A  simplified  theory  of  this  phenomenon, 
which  impairs  the  utilization  of  the  active  material,  was 
presented  by  Dasoyan  and  Aguf  [1],  Several  authors  devel¬ 
oped  for  this  purpose  a  rigorous  mathematical  method  [2-4]. 
We  have  shown  [5]  that  a  cell  of  a  lead-acid  battery  can  be 
modelled  by  an  electrical  equivalent  circuit  consisting  of  a 
pair  of  commercial  lead  grids,  mutually  interconnected  by  a 
system  of  parallel,  thin  resistance  wires.  These  were 
designed  so  as  to  represent  the  sum  of  Ohmic  resistances 
of  electrolyte  and  active  mass  and  of  polarisation  resis¬ 
tances.  The  current  passing  through  each  of  the  wires  was 
determined  from  potential  measurements.  This  approximate 
method  was  used  to  obtain  information  about  the  influence 
of  the  position  of  the  tabs  with  current  leads.  To  obtain 
information  about  the  influence  of  the  charge  passed,  math¬ 
ematical  methods  are  more  suitable.  The  results  are  pre¬ 
sented  here. 
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2.  Equivalent  circuit  and  input  data 

The  mathematical  model  was  based  on  an  electrical 
equivalent  circuit  similar  to  that  used  by  previous  authors 
[3,4].  A  single  mesh  of  this  circuit  is  shown  in  Fig.  1.  It  is 
assumed  that  the  positive  and  negative  grids  have  equal 
numbers  of  meshes  of  equal  forms,  thus  corresponding,  in 
the  geometrical  sense,  to  each  other.  The  nodes  are  visua¬ 
lised  by  black  points.  The  “effective”  Ohmic  resistances  of 
the  horizontal  and  vertical  half  grid  members  including 
active  mass  are  denoted  as  Rx  and  Ry,  respectively,  with 
superscript  (+  or  — )  referring  to  the  plate  sign.  The  internal 
resistance  of  an  elementary  section  of  the  electrode  system 
is  denoted  as  Ry  with  a  whole-numbered  subscript  ( k , 
k  +  1,  etc.).  Similarly,  the  electrode  potentials  at  the  nodes 
are  denoted  as  V  or  W  with  a  whole-numbered  subscript, 
and  the  currents  flowing  through  the  internal  resistances 
RVk,Rvk+l  are  denoted  as  Ik,  4+1.  Thus,  the  electrode  system 
can  be  build  up  by  combining  a  number  of  three-dimen¬ 
sional  meshes  shown  in  Fig.  1.  An  analogous  equivalent 
circuit  was  build  up  from  Ohmic  resistors  in  our  previous 
work  [5]. 

The  resistances,  Ryk,  between  the  nodes,  i.e.  between  the 
electrode  elements  (Fig.  1)  were  now  considered  as  func¬ 
tions  of  the  current,  4,  and  charge  passed,  Qk,  corresponding 
to  Ml  node.  To  obtain  an  analytical  expression  for  this 
function,  the  time  course  of  the  internal  resistance  of  a 
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Nomenclature 

a  length  of  vertical  grid  member 

b  length  of  horizontal  grid  member 

C  discharge  capacity  of  model  cell 
d  thickness  of  active  mass  layer 
Gh  conductance  of  horizontal  grid  member 

G£  conductance  of  horizontal  grid  member  includ¬ 

ing  negative  active  mass 

Gha  conductance  of  horizontal  column  of  negative 
active  mass 

Gv  conductance  of  vertical  grid  member 
Gva  conductance  of  vertical  column  of  negative 
active  mass 

G*  conductance  of  vertical  grid  member  including 
negative  active  mass 
I  total  cell  current 

4  current  flowing  through  internal  resistance  RVk 
Q  total  charge  passed  (It) 

Qk  charge  passed  due  to  current  Ik 

Rx  (with  grid  sign  as  superscript)  Ohmic  resistance 
of  horizontal  half  grid  member 
Ry  (with  grid  sign  as  superscript)  Ohmic  resistance 
of  vertical  half  grid  member 
RVk  Ohmic  resistance  of  kth  elementary  electrode 
system  section 

t  time  of  current  passage 

Uk  voltage  between  Mi  positive  and  negative  nodes 

Vk  electrode  potential  at  Mi  positive  node 

Wk  electrode  potential  at  kth  negative  node 

v  horizontal  coordinate  axis 

y  vertical  coordinate  axis 

o  conductivity  of  negative  active  mass 

laboratory  cell  (including  the  electrolyte,  interphase,  active 
mass,  and  polarisation  resistances)  was  measured  at  various 
discharge  currents.  The  cell  was  assembled  from  five  posi¬ 
tive  and  six  negative  electrodes  of  dimensions  lOmmx 


15  mm  (with  a  Daramic  separator)  cut  out  from  commercial 
SLI  grids  (i.e.  two  meshes  each)  and  connected  in  parallel. 
The  electrolyte  was  H2S04  of  1.28  g/cm3  density.  The  cell 
was  discharged,  in  turn,  with  currents  of  0.1,  0.15,  0.2,  0.3, 
and  0.35  A.  The  internal  cell  resistance  was  determined  at 
10  min.  intervals  by  lowering  the  discharge  current  for  a 
moment  to  one-half,  thus  producing  a  current  drop,  A/,  and 
recording  the  voltage  change,  A  U.  The  sought  resistance 
value  was  set  equal  to  AG/A/.  It  was  taken  into  account  that 
the  resistance  of  the  experimental  cell  was  by  the  factor  of 
0.1  lower  (and  the  current  by  the  factor  of  10  higher)  than 
that  of  the  cell  element  in  Fig.  1.  Thus,  we  arrived  at  the 
value  of  Ryk  corresponding  to  this  cell  element.  The  results 
were  expressed  by  the  following  regression  function  calcu¬ 
lated  by  the  least  squares  method: 

RVk  =  0.32  +  4.2(4  +  2.6  x  10“5  exp(730 Qk  +  1004  -  14) 

(1) 

The  smoothed  dependence  of  RVk  on  Qk  at  various  dis¬ 
charge  currents  is  illustrated  in  Fig.  2.  It  is  apparent  that  the 
most  significant  increase  of  the  resistance  takes  place  when 
the  end  of  discharge  is  approached.  Avery  similar  behaviour 
was  observed  by  Winsel  and  coworkers  [6,7],  who  measured 
the  sum  of  interphase  (contact)  and  active  mass  resistances, 
and  by  Hollenkamp  et  al.  [8],  who  studied  the  interphase 
resistance  between  lead  grid  and  active  mass.  An  exponen¬ 
tial  dependence  of  the  resistance  on  the  degree  of  discharge 
was  also  found  by  Bouet  [9].  It  is  also  apparent  that  the 
resistance  increase  occurs  earlier  at  higher  currents  (Fig.  2). 
This  may  be  attributed  to  the  fact  that  the  current  distribution 
in  the  active  mass  becomes  more  non-uniform  at  higher 
currents,  causing  preferential  discharge  of  the  “outer” 
active  mass  layer.  Accordingly,  the  resistance  of  this  layer 
increases  more  considerably,  and  this  effect  is  reflected  by 
the  behaviour  of  the  Ryk  values. 

With  the  positive  plate,  the  active  mass  does  not  con¬ 
tribute  significantly  to  the  conductance  of  the  grid  members. 
However,  with  the  negative,  the  grid  member  conductances 


Fig.  1.  Part  of  the  equivalent  circuit  used  for  calculation  of  the  current  distribution  over  the  electrode  surface.  Ryk  resistance  between  two  electrode  elements 
(Eq.  (1)).  Positive  grid:  Rf=  1.6125  mQ,  Ry  —  0.5375  mQ;  negative  grid  (with  active  mass  in  the  charged  state):  Rf  =  1.327  mQ,  Rf  =  0.327  mQ;  frame: 
Rx o  =  0.5375  mh,  Ry 0  =  0.3583  mQ. 
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Fig.  2.  Dependences  of  the  internal  resistance  of  a  cell  element  on  the  charge  passed  for  discharge  currents  rising  from  0.01  to  0.035  A.  Experimental  points 
are  fitted  with  regression  curves  based  on  Eq.  (1). 


must  be  corrected  for  the  conductance  of  the  active  mass. 
This  can  be  done  approximately  as  follows:  [3]. 

The  negative  plate  electrode  contains  vertical  grid  mem¬ 
bers  of  length  a  and  horizontal  grid  members  of  length  b,  and 
an  active  mass  layer  of  thickness  d  and  conductivity  a 
(measured  in  Q  1  cm-1).  We  shall  consider  a  column  of 
active  mass  of  thickness  d,  height  a,  and  width  b,  surround¬ 
ing  a  vertical  grid  member  of  length  a.  The  conductance  of 
this  column  in  vertical  direction,  Gva  (measured  in  Q  _1),  is 
given  as 

b 

Gva  =  a  d  -  (2) 

a 

Such  a  column  can  be  assigned  to  each  vertical  grid 
member  except  for  the  vertical  frame  members,  to  which 
only  one-half  of  the  column  of  width  bl 2  belongs.  If  we 
denote  the  conductance  of  the  vertical  grid  member  as  Gv, 
the  conductance  of  this  grid  member  with  the  active  mass 
will  approximately  be  given  by  the  equation 

G*  =  Gv  +  Gva  =  Gy  +  <j  d  -  (3) 

a 

Similarly,  we  may  consider  a  column  of  active  mass  of 
thickness  d,  height  a  and  width  b,  surrounding  a  horizontal 
grid  member  of  length  b.  The  conductance  of  this  column  in 
horizontal  direction,  Gha,  is  given  as 

Gha  =  ad  j  (4) 

b 

One-half  of  such  a  column  of  height  a/2  can  be  assigned  to 
the  horizontal  frame  members.  By  analogy  to  Eq.  (3),  the 
conductance  of  a  horizontal  grid  member  with  the  active 
mass  can  approximately  be  written  as 


The  following  parameter  values  were  used  in  the  calcula¬ 
tions:  a  =  0.5  cm,  b  =  1.5  cm,  d  =  0.1  cm,  a  =  2000  Q  1 
cm-1.  The  latter  value  was  taken  from  the  [3]  and  checked 
by  comparing  the  calculated  potential  drop  within  the 
negative  plate  with  the  measured  one.  The  conductances 
of  the  grid  members  were  determined  as  Gv  =  930  S  and 
Gh  =  310S.  In  accord  with  our  earlier  work  [10],  the 
conductance  of  the  negative  active  mass  decreases  practi¬ 
cally  linearly  with  the  charge  passed,  Q,  during  constant- 
current  discharge;  the  end  value  for  Q  =  0.035  Ah  equals 
0.4  times  the  starting  value.  Thus,  with  the  help  of  the 
definitions  (3)  and  (5),  we  obtain  Gva  =  ad(b/a)  = 
600  —  10  286*2  and  Gha  =  ad{a/b )  =  66.67  —  1142.9 Q. 
Now,  Eqs.  (3)  and  (5)  can  be  put  into  the  form 

G;  =  1530-  10286(9  (6) 

GJ  =  376.67  -  1142.90  (7) 

Hence,  the  dependences  of  the  resistances  in  Fig.  1  on  the 
charge  passed  at  kt h  node,  Qk,  can  be  expressed  by  the 
following  approximate  equations: 


R-  _  (  °-5 
*  V376.67  -  1142.90* 

_  _  /  0.5 

Ry  ~  V1530^1o286^ 

where  Qk  =  f  /*  dr  and  /*  is  generally  different  for  different 
nodes;  in  other  words,  these  two  quantities  depend  on  the  jc 
and  y  coordinates  (used  in  3D  diagrams,  Figs.  3-8). 


3.  Method  of  calculation 


G\ i  +  Gha  =  Gh  +  a  d  - 

b 


(5) 


By  applying  the  first  and  second  Kirchhoff  laws  to  the 
nodes  and  loops  of  the  whole  equivalent  circuit,  a  system  of 
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Fig.  3.  Current  distribution  over  the  plate  electrode  surface  at  increasing  state  of  discharge,  Q,  from  0  to  100%  for  the  most  unfavourable  tab  configuration 
(units  on  x  and  y  axes  are  arbitrary  in  this  and  subsequent  figures). 


linear  equations  is  obtained,  whose  solution  gives  the  sought 
distribution  of  potentials  and  currents.  To  take  into  account 
the  system  changes  with  the  time,  the  calculations  were 
carried  out  in  steps  as  follows. 


the  internal  resistance  RVk  and  the  potential  distribution  in 
the  nodes  of  the  equivalent  circuit  is  calculated.  The  node 
potentials  of  Ath  element,  V]k  and  W]k ,  are  used  to  calculate 
the  corresponding  voltage  and  current 


•  Step  1  (/  =  1,  t\  —  1  s,  Ati  =  1  s):  The  initial  value  of 
R[  =  0.32  Q  is  introduced  into  all  elements  representing 


(10) 
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Fig.  4.  Current  distribution  over  the  plate  electrode  surface  at  increasing  state  of  discharge,  Q,  from  0  to  100%  for  the  tab  configuration  shown. 


The  charge  passed  through  the  Mi  element  can  be  calcu¬ 
lated  from  the  current  as 

(ID 

Step  2  (/  =  2,  t2  =  30  s,  f 2  ~  0):  The  internal 

resistance  corresponding  to  the  Mi  element  is  calcu¬ 
lated  from  Eq.  (1).  Afterwards,  the  node  potential 


distribution  is  again  calculated  and  the  corresponding 
voltage  and  current  values  are  found  as  in  step  1.  The 
charge  passed  through  the  Mi  element  is  then  calculated 
as 

Ql  =  Ql+i2k/Sh  (12) 

The  procedure  continues  step  by  step  until  tn  =  12,000  s. 


0  =  60  % 


/,  mA  24' 


0  =  80  % 


/,  mA  2 


0  =  90  % 

27- 

/,  mA  24‘ 


zSziSSSgsi 

11^  1111111 


Wsj0&£ 


0=100 


/,  mA  2 


C^I 


X 

Fig.  6.  Current  distribution  over  the  plate  electrode  surface  at  increasing  state  of  discharge,  2,  from  0  to  100%  for  the  tab  configuration  shown. 


and  100%  discharge.  This  corresponds  to  discharge  times  of 
0,  80,  120,  160,  180,  and  200  min,  discharge  current 
/  =  2  A,  and  discharge  capacity  C  =  6.66  Ah.  Six  config¬ 
urations  of  the  current  tabs  welded  to  the  lead  grids  were 
considered.  These  are  indicated  schematically  in  the  dia¬ 
grams. 


It  can  be  seen  that  the  initial  rate  of  discharge  (i.e.  the  local 
current)  is  always  more  or  less  elevated  in  regions  close  to 
the  tabs.  This  phenomenon  is  most  pronounced  in  Fig.  3, 
where  the  tabs  are  located  in  the  corners  facing  each  other. 
Therefore,  this  region  reaches  the  discharged  state  more 
rapidly.  Since  the  total  current  is  maintained  constant,  the 
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Fig.  7.  Current  distribution  over  the  plate  electrode  surface  at  increasing  state  of  discharge,  Q,  from  0  to  100%  for  the  tab  configuration  shown. 


current  drop  in  one  region  causes  the  discharge  rate  in 
another  one  to  increase,  although  not  as  much,  as  can  be 
seen  from  the  diagrams.  In  spite  of  this  levelling  effect,  the 
situation  at  the  end  of  discharge  for  most  tab  configurations 
(Figs.  3-7)  does  not  look  satisfactory.  A  moderate  improve¬ 
ment  can  be  seen  from  Fig.  7,  but  the  best  configuration 


seems  to  be  the  last  one,  with  extended  tabs  on  opposite 
ends  of  the  plates,  corresponding  to  Fig.  8.  This  is  in  accord 
with  the  theory  presented  in  our  preceding  work  [5].  Here, 
the  current  distribution  over  the  plate  surface  remains 
fairly  uniform  all  the  time,  with  a  shallow  extremum  in 
the  middle  region  of  the  plates.  Since  the  conductivity  of 
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Fig.  8.  Current  distribution  over  the  plate  electrode  surface  at  increasing  state  of  discharge,  Q,  from  0  to  100%  for  the  most  favourable  tab  configuration. 


the  negative  active  mass  was  taken  into  account,  the  negative 
plate  is  a  better  conductor  than  the  positive,  resulting  in 
a  shift  of  the  position  of  the  extremum  toward  the  tab  of 
the  negative  electrode  (compare  the  theory  in  our  previous 
work  [5]). 

An  illustrative  criterion  characterising  the  uniformity  of 
the  current  distribution  is  the  ratio  of  the  maximum  to  the 
minimum  local  currents  (ImdX/Im\n),  which  is  plotted  as  a 
function  of  the  depth  of  discharge,  Q,  in  Fig.  9  for  all  the  tab 


configurations  under  discussion.  The  configuration  used  by 
many  manufacturers  takes  the  mean  position.  The  config¬ 
uration  with  extended  tabs  located  on  opposite  ends  of  the 
positive  and  negative  plates  is  obviously  by  far  the  best  one, 
since  the  mentioned  criterion  is  only  a  little  higher  than  one. 
Recently,  Lam  and  coworkers  [11]  proposed  a  somewhat 
similar  but  more  complicated  arrangement,  where  the  tabs 
are  located  at  the  top  and  bottom  of  each  of  the  positive  and 
negative  plate  electrodes. 
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Fig.  9.  Dependences  of  the  criterion  /max//mjn  on  the  relative  state  of  discharge  for  seven  configurations  of  tabs  welded  to  the  lead  grids. 
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